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ABSTRACT 

We present results of deep echelle spectrophotometry of the brightest knot of the 
Herbig-Haro object HH 202 in the Orion Nebula -HH 202-S- using the Ultraviolet 
Visual Echelle Spectrograph (UVES) in the spectral range from 3100 to 10400 A. The 
high spectral resolution of the observations has permitted to separate the component 
associated with the ambient gas from that associated with the gas flow. We derive 
electron densities and temperatures from different diagnostics for both components, 
as well as the chemical abundances of several ions and elements from collisionally ex- 
cited lines, including the first determinations of Ca + and Cr + abundances in the Orion 
Nebula. We also calculate the He + , C 2+ , + and 2+ abundances from recombina- 
tion lines. The difference between the 2+ abundances determined from collisionally 
excited and recombination lines -the so-called abundance discrepancy factor- is 0.35 
dex and 0.11 dex for the shock and nebular components, respectively. Assuming that 
the abundance discrepancy is produced by spatial variations in the electron tempera- 
ture, we derive values of the temperature fluctuation parameter, t 2 , of 0.050 and 0.016, 
for the shock and nebular components, respectively. Interestingly, we obtain almost 
coincident t 2 values for both components from the analysis of the intensity ratios of 
He I lines. We find significant departures from case B predictions in the Balmer and 
Paschen flux ratios of lines of high principal quantum number n. We analize the ion- 
ization structure of HH 202-S, finding enough evidence to conclude that the flow of 
HH 202-S has compressed the ambient gas inside the nebula trapping the ionization 
front. We measure a strong increase of the total abundances of nickel and iron in the 
shock component, the abundance pattern and the results of photoionization models for 
both components are consistent with the partial destruction of dust after the passage 
of the shock wave in HH 202-S. 

Key words: ISM: abundances - Herbig-Haro objects - ISM: individual: Orion Nebula 
- ISM: individual: HH 202 



1 INTRODUCTION 

HH 202 is one of the brightest and most conspicuous Herbig- 
Ha ro (HH) obje c ts of the Orion Nebula. It was discovered 
by ICanto et all (|l980T ). The origin of this outflow is not 
clear, though the radial velocity and proper motion stud- 
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ies suggest that this object forms a great complex together 
with HH 203, 204, 269, 529, 528 and 625, with a com- 
mon origin in one or more sources embedd ed within the 
Orion Molecular Cloud 1 South (OMC IS) (seelRosado et ail 



2002: IO'Dell fc Doi 



20031 : lO'Dell fc Henn cv 2008). Recently, 



iHennev et al.l (|2007j ) have summarized the main character- 
istics of these outflows a nd an extensive study of their kine- 
matics can be found in iGarcia-Diaz et al.l (2008). HH 202 
shows a wide parabolic form with several bright knots of 
which HH 202-S is the brightest one (see Figure [T]). 
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Figure 1. HST image of the central part of the Orion Nebula combined with WFPC2 images in different filters jO'Dell &: Wong 1996). 
The white square corresponds to the field of view of the Potsdam Multi-Aperture Spectrograph (PMAS) and the separate close up image 
at the right shows the rebinned Ha map presented in Paper I. Inside of this box, the black rectangle indicates the slit position and the 
area covered by the UVES spectrum analyzed in this paper (l'.'5x2'.'5). 



The kinematic properties of HH 202-S have been studied 
by means of high-spectral re solution spectroscopy by sev- 
eral authors. iDoi et all (|2004l ) have found a radial velocity 
of — 39±2 km s ~ 1 , i n agreement with p revious results by 
Mca burnl l|l986l ) and lO'Dell et al.l (|l99lh . lO'Dell fc Hennevl 
1 2008) have determined a tangential velocity of 59±8 
k m s -1 , which is i n agre e ment with previou s deter minations 
by lO'Dell fc Doil (|2003l ). lO'Dell fc Hennevl |2008l ) have cal- 
culated a spatial velocity of 89 km s _1 and an angle of the 
velocity vector of 48° with r espect to the plane o f the sky, 
similar to the values found bylHennev et al] (|2007l ). Imaging 
studies bv lO'Dell et~ai1 1| 19971 ) with the Hubble Space Tele- 
scope (HST) of the HH objects in the Orion Nebula show 
an extended [O ill] emission in HH 202 and strong [O ill] 
emission in HH 202-S. This fact, together with the closeness 
of HH 202 to the main ionization source of the Orion Neb- 
ula, 9 1 Ori C, indicate that the excitation of the ionized gas 
is dominated by photoionization in HH 202-S, though the 
observed radial velociti es imply that som e shocked gas can 
be mixed in the region (|Cant6 et al.lll980l ). Photoionization- 
dominated flows are a minority in the inventory of HH ob- 
jects, which are typically excited by shocks. This kind of HH 
objec ts is also known as "irradiated jets" (|Reipurth et al.l 
1998), since they are excited by the UV radiation from 
nearby massive st ars. Irradiated jets have been found in the 
Orion Nebula (e.g-lBallv fc Reipurthl200ll;lBallv et al.ll2006l; 
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20031). the Carina Nebula (|Smith et all .20041) . NGC 1333~ 
(Ball y et alj|2006l) and the Trifid Nebula ([Cernicharo et~ail 



19981: iReipurth et al l 



1998). 



Mesa-Delgado et al.l (j2008a|) have obtained the spatial dis- 



tributions of the physical conditions and the ionic abun- 
dances in the Orion Nebula using long-slit spectroscopy 
at spatial scales of 1'.'2. The goal of that work was to 
study the possible correlations between the local struc- 
tures observed in the Orion Nebula -HH objects, pro- 
ply ds, ionization fronts- and the abundance discrepancy 
(AD) that is found in H II regions. The AD is a clas- 
sical problem in the study of ionized nebulae: the abun- 



dances of a given ion derived from recombination lines, 
RLs, are often between 0.1 and 0.3 dex higher than those 
obtained fro m collisionally excited lines, CELs, in H H 
regio n s fsee iGarcfa-Rojas fc Estebanl 120071 : lEsteban et alJ 
120041 : iTsamis et all l2003h . The difference between those 
independent determinations of the abundance defines the 
abundance discrepancy factor, ADF. The predicti ons of the 
tempe rature fluctuation paradigm proposed by IPeimbertl 
l|l967t . -and parametrized by the mean square of the spa- 
tial distribution of temperature, the t 2 parameter- seem to 
account for the discrepancies observed in H 11 regions (see 
iGarcfa-Roias fc Estebanl l2007t ) . A striking result found in 
the spatially-resolved study of iMesa-Deleado et al.l (|2008al ) 
is that the ADF of 2+ , ADF(0 2+ ), shows larger values at 
the locations of HH objects as is the case of HH 202. Using 
integral field spectroscopy with interme diate-spectral resolu- 
tion an d a spatial resolution of l"xl", Me sa-Delgado et al.l 
(2008b|) (hereinafter, Paper I) have mapped the emission 
line fluxes, the physical properties and the 2+ abun- 
dances derived from RLs and CELs of HH 202. They have 
found extended [O ill] emission and higher values of the 
electron density and temperaure as well as an enhanced 
ADF(Q 2+ ) in HH 202-S, c onfirming the earlier results of 
iMesa-Deleado et al.1 l|2008ah . 

HH 529 is another HH object that is photoionized by 8 1 
Ori C and shows simi lar characteristics to those of HH 202. 
iBlagrave et al] (|2006l ) have performed deep optical echelle 
spectroscopy of that object with a 4m-class telescope and 
have detected and measured about 280 emission lines. Their 
high-spectral resolution spectroscopy allowed them to sepa- 
rate the kinematic components associated with the ambient 
gas and with the flow. They have determined the physical 
conditions and the ionic abundances of oxygen from CELs 
and RLs in both components. However, they do not find 
high ADF(0 2+ ) and t 2 values in neithe r comp onent. An- 
other interesting result of IBlagrave et all (|2006l ) is that the 
ionization structure of HH 529 indicates that it is a matter- 
bounded shock. 

Motivated by the results found by iMesa-Deleado et al.l 
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Figure 2. Sections of the bidimensional UVES spectrum showing 
the spatio-kinematic profiles of different oxygen lines: [O i] 6300 
A, O I 7772 A, [O II] 3728 A, O n 4649 A and [O in] 4959 A. 
Wavelength increases to the right and north points up. The solid 
straight lines in all the sections correspond to the slit center and 
the dashed lines —only represented in the case of the [O II] 3728 
A profile— correspond to the extracted area 2'.'5 wide. 



l|2008al ). inspired by the work of iBlagrave et all l|2006l ) and 
in order to complement the results presented in Paper I, 
we have isolated the emission of the flow of HH 202-S knot 
using high-spectral resolution spectroscopy, presenting the 
first complete physical and chemical analysis of this knot. 
In [J2] we describe the observations of HH 202 and the re- 
duction procedure. In jj3]we describe the emission line mea- 
surements, identifications and the reddening correction, we 
also compare our reddening determinations with those avail- 
able in the literature. In £0]we describe the determinations 
of the physical conditions, the chemical -ionic and total- 
abundances and the ADF for + and 2+ . In Sj5]we discuss: 
a) some inconsistencies found in the Balmer decrement of 
the lines of higher principal quantum number, b) the ion- 
ization structure of HH 202-S, c) the radial velocity pattern 
of the lines of each kinematic component, d) the t 2 parame- 
ter obtained from different methods and its possible relation 
with the ADF, and e) the evidences of dust grain destruc- 
tion in HH 202-S. Finally, in fJS] we summarize our main 
conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

HH 202 was observed on 2003 March 30 at Cerro Paranal 
Observatory (Chile), using the UT2 (Kueyen) of the Very 
Large Telescope (VL T) with the Ultraviolet Visual Echelle 
Spectrograph fUVES. lLvbdorico et aI1l2000h . The standard 
settings of UVES were used covering the spectral range from 
3100 to 10400 A. Some narrow spectral ranges could not 



be observed. These are: 5783-5830 and 8540-8650 A, due 
to the physical separation between the CCDs of the de- 
tector system of the red arm; and 10084-10088 and 10252- 
10259 A, because the last two orders of the spectrum do 
not fit completely within the size of the CCD. Five individ- 
ual exposures of 90 seconds -for the 3100-3900 and 4750- 
6800 A ranges- and 270 seconds -for the 3800-5000 and 
6700-10400 A ranges- were added to obtain the final spec- 
tra. In addition, exposures of 5 and 10 seconds were taken 
to obtain good flux measurements -i.e. non-saturated- for 
the brightest emission lines. The spectral resolution was 
A/AA ~ 30000. This high spectral resolution enables us 
to separate two kinematic components: one corresponding 
to the ambient gas -which we will call nebular component 
and whose emission mainly arises from behind HH 202 and, 
therefore, could not entirely correspond to the pre-shock 
gas- and another one corresponding to the gas flow of the 
HH object, the post-shock gas, which we will call shock com- 
ponent. 

The slit was oriented north-south and the atmospheric dis- 
persion corrector (ADC) was used to keep the same observed 
region within the slit regardless of the air mass value. The 
HH object was observed between airmass values of 1.20 and 
1.35. The average seeing during the observation was f/.'7. The 
slit width was set to l'.'5 as a compromise between the spec- 
tral resolution needed and the desired signal-to-noise ratio 
of the spectra. The slit length was fixed to 10". The one- 
dimensional spectra were extracted for an area of l'.'5x2 / .'5. 
This area covers the apex of HH 202, the so-called knot 
HH 202-S, as we can see in Figure [T] This zone shows the 
maximum shift in velocity between the shock and nebular 
components (see Figure [2} allowing us to appropriately sep- 
arate and study the spectra of both kinematic components. 
The spectra were reduced using the iRAfQ echelle reduction 
package, following the standard procedure of bias subtrac- 
tion, aperture extraction, flatfielding, wavelength calibration 
and flux calibration The standard stars EG 247, C-32 9927 
ilHarmiv et alll99l[l994l) a nd HD 49798 (|Bohlin fc Lindlerl 
1993: iTurnshek et alj Il990h were observed to perform the 
flux calibration. The error of the absolute flux calibration 
was of the order of 3%. 



3 LINE MEASUREMENTS, 

IDENTIFICATIONS AND REDDENING 
CORRECTION 

Line fluxes were measured applying a double Gaussian pro- 
file fit procedure over the local continuum. All these mea- 
surements were made with the SPLOT routine of iraf. 
All line fluxes of a given spectrum have been normalized 
to a particular bright emission line present in the common 
range of two consecutive spectra. For the bluest spectrum 
(3100-3900 A), the reference line was H9 3835 A. For the 
range from 3800 to 5000 A, the reference line was H/3. In 
the case of the spectrum covering 4750-6800 A, the refer- 
ence was [O in] 4959 A. Finally, for the reddest spectrum 
(6700-10400 A), the reference line was [S n] 6731 A. In or- 
der to produce a final homogeneous set of line flux ratios, 

1 IRAF is distributed by NOAO, which is operated by AURA, 
under cooperative agreement with NSF. 
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Figure 3. Section of the echelle spectrum of HH 202-S showing 
the shock (left) and nebular (right) components of each emission 
line of multiplet 1 of O II. 



all of them were rescaled to the H/3 flux. In the case of the 
bluest spectra the ratios were rescaled by the H9/H/3 ratio 
obtained from the 3800-5000 A ranee. The emission line ra- 
tios of the 4750-6800 A range were multipled by the [O ill] 
4959/H/3 ratio measured in the 3800-5000 A range. In the 
case of the last spectral section, 6700-10400 A, the rescal- 
ing factor was the [S n] 6731/H/3 ratio obtained from the 
4750-6800 A spectrum. All rescaling factors were measured 
in the short exposure spectra in order to avoid the possible 
saturation of the brightest emission lines. This process was 
done separately for both the nebular and shock components. 
The spectral ranges present overlapping regions at the edges. 
The adopted flux of a line in the overlapping region was ob- 
tained as the average of the values obtained in both spec- 
tra. A similar procedure was considered in the case of lines 
present in two consecutive spectral orders of the same spec- 
tral range. The average of both measurements was consid- 
ered for the adopted value of the line flux. In all cases, the 
differences in the line flux measured for the same line in dif- 
ferent orders and/or spectral ranges do not show systematic 
trends and are always within the uncertainties. 
The identification and laboratory wavelengths of the lines 
we re obtained foll o wing a previous work on the Orion Nebula 
bv lEsteban et al.l (|2004l ). the compilations by iMoorel (|l945l ) 
and The Atomic Line List v2.o4E The identification process 
and the measurement of line fluxes were done simultane- 
ously. The inspection of the line shapes at the bi-dimensional 
echelle spectrum was always used to identify which compo- 
nent -nebular or shock- was measured at each moment. The 
rather different spatial and spatio-kinematic structure of the 
two kinematic components is illustrated in Figure[2] We have 
identified 360 emission lines in the spectrum of HH 202-S, 
115 of them only show one component -8 belong to the neb- 
ular component and 107 belong to the shock one- and 8 are 
dubious identifications. 

For a given line, the observed wavelength is determined by 
the centroid of the Gaussian fit to the line profile. For lines 
measured in different orders and/or spectral ranges, the av- 
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erage of the different wavelength determinations has been 
adopted. From the adopted wavelength, the heliocentric ve- 
locity, Vhei, has been calculated using the heliocentric cor- 
rection appropriate for the coordinates of the object and the 
moment of observation. The typical error in the heliocentric 
velocity measured is about 1-2 km s _1 . 
All line fluxes with respect to H/3, F(\)/F(H.f3), were dered- 
dened using the typical relation, 



J(A) 

/(H/3) F(H0) 



F W 10 c(H/9)/(A) 



(1) 



The reddening coefficient, c(H/3), respresents the amount of 
interstellar extinction which is the logarithmic extinction at 
H/3, while /(A) is the adopted extinction curve normalized 
to /(H/3) = 0. The reddening coefficient was determined 
from the comparison of the observed flux ratio of Balmer 
and Paschen lines -those not contaminated by telluric or 
other nebular emissions- with respect to H/3 and t he theo- 
retical ones computed bv lStorev fe Hummer! l|l995h for the 
physical conditions of T e = 10000 K and n e = 1000 cm -3 . 
As in Paper I, we have used the reddening functio n, /(A), 
normalized to H/3 derived bv lBlagrave et all (|2007T l for the 
Orion Nebula. The use of this extinction law instead of the 
classical one (Cos tero fc Peimberdfl970l ) produces slightly 
higher c(H/3) values and also slightly different dereddened 
line fluxes depending on the spectral range (see Paper I). 
The final c(H/3) values obtained for the two kinematic com- 
ponents were weighted averages of the values obtained for 
the individual lines: c(H/3) ncb = 0.41±0.02 and c(H/3) sh = 
0.45±0.02. Although not all the c(H/3) values are consistent 
with each other (see H5.1J1 . the average values obtained are 
quite similar and consistent within the uncertainties. 
We can compare the reddening values with those obtained 
from integral field spectroscopy data presented in Paper I in 
the same area of HH 202-S (see Figure [TJ and corresponding 
to the section Aa = [—4, —6] and AS = [—8, —5] (see Figure 
3 of Paper I). The average c(H/3) in this zone is 0.65±0.15, 
which is higher than those determined for UVES data. How- 
ever, if we re-calculate the value of c(H/3) from the UVES 
data using the same Balmer lines as in Paper I, we obtain a 
value 0.5±0.1 in both kinematic components, a value consis- 
tent with the PMAS one within the errors. These differences 
can be related to several systematical disagreements found 
between the c(H/3) values obtained from different individual 
Balmer or Paschen lines (see H5.ip . 

In the most complete w ork on the reddening d i strib ution 
across the Orion Nebula, lO'Dell fc Yusef-Zadehl (|200Ch ob- 
tain values of c(H/3) between 0.2 and 0.4 in the zone around 
HH 202-S, somewhat lower than our reddening determina- 
tions. This can be due to the fact that O'Dell fc Yusef-Zadehl 
use the extinction law bv lCostero fc Peimbertl <|l970l ) which, 
as we discuss in Paper I, produc es lower c(H/3) values than 
the more recent extinction law dBlagrave et al.l|2007l ). We 
have also re- calculated c(H/3) from our UVES spectra mak- 
ing use of the lCostero fc Peimbertl law, and we obtain values 
a bout 0.3, being now in agreem ent with the determinations 
of lO'Dell fc Yusef-Zadehl (|2000h . 

In Table [TJ the final list of line identifications (columns 1— 
3), /(A) values (column 4), heliocentric velocities (columns 
5 and 8) and dereddened flux line ratios (columns 6 and 9) 
for the nebular and shock component are presented. The ob- 
servational errors associated with the line dereddened fluxes 
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with respect to H/3 -in percentage- are also presented in 
columns (7) and (10) of Table [1] These errors include the 
uncertainties in line flux measurement, flux calibration and 
error propagation in the reddening coefficient. 
In column (11) of Table [1] we present the shock-to-nebular 
line flux ratio for those lines in which both kinematic com- 
ponents have been measured. This ratio is defined as: 

7 sh = [7(A)/7(H/3)] sh 

-^ncb 



[7(A)/7(H/3)] n 



7(A) sh 7(H/3) ncb 



7(A) ne 



7(H/3) S 



(2) 



where the integrated dereddened H/3 fluxes are 7(H/3) nc b 
= (3.80±0.20)xl0~ 12 erg cm" 2 s" 1 and 7(H/3) sh = 
(6.00±0.20)xl0" 12 erg cm" 2 s" 1 . The 7 ncb /7 sh ratios de- 
pend on each particular line. In general, they are close to 1 
for H I lines but become less than 1 for higher ionized species 
-except Fe ions-, and are typically greater than 1 for neutral 
species. A more extensive discussion on this particular issue 
will be presented in section t|5.2l 

In Figure [3] we show a section of our flux-calibrated echelle 
spectra around the lines of multiplet 1 of O II. It can be 
seen that both the nebular and the shock components are 
well separated and show a remarkable high signal-to-noise 
ratio. 



4 RESULTS 

4.1 Excitation mechanism of the ionized gas in 
HH 202-S 



In their recent work, lO'Dell fe Hennevl l|2008l ) argue that 
the presence of a variety of ionization stages in the ionized 
gas of HH 202 indicates that the flow also contains neutral 
material. They interpret that fact as due to the impact of 
the flow with preexisting neutral material -perhaps of the 
foreground veil- or that the flow compresses the ambient 
ionized gas inside the nebula to such degree that it traps 
the ionization front. Our results provide some clues that 
can help to ascertain this issue. The value of some emission 
line ratios are good indicators of the presence of shock ex- 
citation in ionized gas, especially [S ll]/Ha and [O i]/Hq. 
In our spectra, we find log([S n] 6717+31/Ha) values which 
are almost identical in both kinematic components ( — 1.49 
and —1.44 for the nebular and shock component, respec- 
tively). These values are completely consistent with those 
expected for photoionized nebulae and far from the range of 
values between —0.5 and 0.5, which is the typic al of super- 
nova remnants and HH objects (see Figure 10 o f R iera et al.l 
1 19891 ). On the other hand, the values of log([0 i] 6300/Ha) 
that we obtain for the nebular and shock component are of 
—2.66 and —2.22, somewhat different in this case, but also 
far from the values expected in the case of substa ntial con- 
tribution of shock excitation (|Hartigan et al.|[l987 | ). Finally, 
we ha ve also used the diagnostic diagrams of Raga et al.l 
|2008l ) where the [N n] 6548/Hq and [S II] 6717+31/Hq vs. 
[O III] 5007/Hq ratios of HH 202 are found in the zone dom- 
inated by photoionized shocks. Therefore, the spectrum of 
HH 202-S seems to be consistent with the picture that the 
bulk of the emission in this area is produced by photoion- 
ization acting on compressed ambient gas that has trapped 
the ionization front inside the ionized bubble of the nebula. 
In the rest of the paper we will provide and discuss further 
indications that HH 202-S contains an ionization front. 



4.2 Physical Conditions 

We have computed physical conditions of the two kine- 
matic components using several ratios o f CELs following the 
same m ethodology as in Paper I and in Mesa-Delgado et al.l 
(2008a]). The electron temperatures, T e , and densities, n e , 
are presented in Tabled 

We have determined n e from [O n], [S n], [ CI hi] and [Ar iv] 
line ratios using the nebular package (Sha w fe Dufourl 
Il995h . In the case of the n e obtained from [Fe ill] lines, 
we have used flux ratios of 31 and 12 lines for the shock 
and nebular comp onent, respectively, follow ing the proce- 
dure described by iGarcfa- Roias et all l|2006h . For the neb- 
ular component, we have adopted the average value of 
n e ([0 n]), n e ([S n]) and n e ([Cl ill]) excluding n e ([Fe ill]) 
and n e ([Ar iv[) due to their discrepant values and very large 
uncertainties. For the shock component, we have adopted 
the average of n e ([0 n]), n e ([Cl ill]) and n e ([Fe ill]), while 
n e ([S ll[) has not been included because the [S n] line ratio 
is out of the range of validity of the indicator. As we can 
see in Table [2] the density of the shock component (~17000 
cm" 3 ) is much higher than the density of the nebular one 
(-3000 cm" 3 ). 

However, the bulk of the emission of the nebular component 
might come from behind HH 202, and the electron density 
that we have found for that component might not be the true 
one of the pre-shock gas. In fact, tak ing into account that the 
velocity of the gas flow is 89 km s" 1 (|Q'Dell fe Hennevll2008f ) 
and the typical sound speed of an ionized gas is about 10-20 
km s , we have adopted a Mach number, M, for HH 202 
of about 5 and, thus, the shock compression ratio should be 
M 2 ~ 25. Using the density of the shock component (see 
Table 0, we obtain a pre-shock density ~ 17430/25 « 700 
cm" 3 . This value is lower than the 2890 cm" 3 determined 
for the nebular component. Therefore, it seems clear that 
the bulk of the nebular component does not refer to the gas 
in the immediate vicinity of HH 202 as we have mentioned 
in Sj2] 

Electron temperatures have been derived from the classi- 
cal CEL ratios of [N n], [O n], [S n], [O III], [S III] and 
[Ar in]. Under the two-zone ionization scheme we have 
adopted T e ([N n]) as representative for the low ionization 
zone and T e ([0 Hi]) for the high ionization zone. We have 
also d erived T e (He i) using the method of iPeimbert et all 
(|2002f ) and state-of-the-art atomic data (see M4.4p . 
We have compared these temperature determinations with 
those obtained from the integral field unit (IFU) data pre- 
sented in Paper I. We have determined the mean T e values 
of the spaxels of the section of the FOV of the PMAS data 
that encompasses the area covered by our UVES spectrum, 
finding <T e ([0 m])> = 8760±260 K and <T C ([N n])> = 
9730±590 K. These values are in agreement within the er- 
rors with those obtained in this paper (see Table |5). The av- 
erage density from the IFU data, obtained from the [S n] line 
ratio, is 7300±3000 cm -3 , a value between the n e adopted 
for each kinematic component from the UVES data. 
As we can see in Table [5J the T e values are quite similar 
in both components with differences of the order of a few 
100 K. The temperatures derived from [N n] lines are higher 
than those derived from [O ill] lines, which is a typical re- 
sult observed in previous works on the Orion Nebula (e.g. 
iMesa-Delgado et al.ll2008al ; iRubin et al.ll2003h . as well as in 
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Table 1. Identifications, reddening-corrected line ratios (I(H/3)=100) for an area of l'.'5x2'/5 and heliocentric velocities 
for the nebular and shock components. 



Nebular Component Shock Component 



Ao (A) a 
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Mult a 


/(A) 6 




/(A) d 


Error (%) e 




I(X) d 


Error (%) e 
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[Fe in] 
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0. 
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0. 
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0. 
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[Fe III] 
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0. 
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0. 
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Table 1. continued 
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Table 1. continued 



Nebular Component Shock Component 

A (A) a Ion" Mult" /(A) 6 V he f I(X) d Error (%) e V hd c I(X) d Error (%) e I ah /I ne b f Notes 
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Nebular Component Shock Component 
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Table 1. continued 



Nebular Component Shock Component 
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Table 1. continued 
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-34 


0.019 


18 






9903.46 


C 11 


17.02 


-0.557 


14 


0.048 


40 


-33 


0.045 


40 


0.947 




9960.00 


[Fc in] 


8F 


-0.562 








-42 


0.022 


30 






10027.7 


He 1 


6/7 


-0.567 


14 


0.152 


9 


-36 


0.159 


9 


1.048 
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Table 1. continued 



Nebular Component Shock Component 



Ao (A) a 


Ion° 


Mult a 




v hel c 


I(X) d 


Error (%) e 


v hel c 




Error (%) e 


hh/Incb f Notes 


10031.2 


He I 


7/7 


-0.568 


11 


0.045 


15 


-38 


0.054 


11 


1.185 


10049.4 


H I 


P7 


-0.569 


13 


5.425 


7 


-38 


5.016 


7 


0.924 


10320.5 


[S II] 


3F 


-0.590 


22 


0.194 


9 


-37 


1.016 


7 


5.234 


10336.4 


[S ii] 


3F 


-0.591 


22 


0.216 


9 


-38 


0.891 


7 


4.117 


10370.5 


[S ii] 


3F 


-0.593 








-38 


0.396 


9 





a Identification of each line: laboratory wavelength, ion and multiplct. 

Value of the extinction curve adopted l lBlagrave et alj|2007l) . 
c Heliocentric velocity in units of km s _1 , the typical error is 1-2 km s _1 . 
d Dereddened fluxes with respect to /(H/3) = 100. 

e Error of the dereddened flux ratios. Colons indicate errors larger than 40 per cent. 

f Shock-to- nebular line flux ratio. See definition in equation (2). 

9 Line blended with another line and deblended via Gaussian fitting. 

h Contaminated by "ghost". 

1 Contaminated by telluric emissions and not deblended. 
J Deblended from telluric emissions. 



Table 2. Physical Conditions. 







Nebular 


Shock 


Indicator 


Component 


Component 


n e (cm 3 ) 


[On] 


3490 ± 810 


18810 ± 8280 




[S n] 


2350 ± 910 


>14200 




[CI in] 


2470 ± 1240 


23780 ± 13960 




[Fe III] 


11800 ± 9000 


17100 ± 2500 




[Ar iv] 


5800 : 






adopted 


2890 ± 550 


17430 ± 2360 


Te (K) 


[N II] 


9610 ± 390 


9240 ± 300 




[O ii] 


8790 ± 250 


9250 ± 280 




[S ii] 


8010 ± 440 


8250 ± 540 




[O in] 


8180 ± 200 


8770 ± 240 




[S in] 


8890 ± 270 


9280 ± 300 




[Ar in] 


7920 ± 450 


8260 ± 410 




He 1 


8050 ± 150 


7950 ± 200 



Paper I. This is a likely result of the ionization stratification 
in the nebula. It is interesting to note that the difference be- 
tween both temperatures is smaller in the case of the shock 
component, in this case, all the emission comes from a - 
probably- much narrower slab of ionized gas. 
The relatively low uncertainties in the physical conditions 
are due to the high signa l-to-noise ratio of the emission lines 
used in the diagnostics. iBlagrave et al.l (|2006l ) computed the 
physical conditions for HH 529 and they obtained similar re- 
sults -higher densities in the shock component but similar 
temperatures in both components- though with compara- 
tively larger errors. 

4.3 Ionic abundances from CELs 

Ionic abundances of N+, 0+, 2+ , Ne 2+ , S+, S 2+ , C1+, 
Cl 2+ , Ar 2+ and Ar 3+ have been derived from CELs un- 
der the two-zone scheme and t 2 = 0, using the nebular. 
package. All abundances were calculated for the shock and 
nebular component, except for Ar 3+ , which was not de- 
tected in the spectrum of the shock component. The atomic 



data for G1+ are not implemented in the nebular routines, 
so we h ave used an old versio n of the 5-level atom pro- 
gram of IShaw fe Dufourl d 19951 ) -FIVEL- that is described 
bv lDe Robertis et all (|l987T ). This prog r am u ses the atomic 
data for this ion compiled by iMendozal (119831 ) . 
We have also measured [Ca n], [Cr n], [Fe Ii], [Fe III], [Fe iv], 
[Ni ii] and [Ni ill] lines. The abundances of these ions are 
also presented in Table They were computed assuming 
the appropriate temperature under the two- zone scheme and 
the procedures indicated below. In addition and only in the 
shock component, we have detected a substantial number 
of lines of other quite rare heavy-element ions as [Cr III], 
[Co ii], [Co III], [Ti ill] and, possibly, [Cr iv], [Co rv], [Mn n] 
and [V ii] . Unfortunately, we cannot derive abundances from 
these lines due to the lack of atomic data for these ions. 
Two [Ca ii] lines at 7291 and 7324 A were detected in the 
shock component. In order to derive the Ca + abundance, we 
solved a 5-level model atom using the singl e atomic dataset 
available for this ion |Melendez et al.ll2007l ). Note that this 
is the first determination of the Ca + abundance in the Orion 
Nebula and this poses a lower limit to the gas-phase Ca/H 
ratio in this object. 

Two and four [Cr n] lines were measured in the nebular and 
shock component, respectively, although those at 8309 and 
8368 A are very faint. [Cr n] lines can be affected by con- 
tinuum or starlight fluorescence as is also the case for the 
[Fe ii] and [Ni n] lines. We have computed the Cr + abun- 
dances using a 180-level mo del atom that t r eat co ntinuum 
fluorescence excitation as in Bautista et al.l (1 19961 ) and in- 



cludes the atomic data of iBautista et all ( 20081 ). In order to 



consider the continuum fluorescence excitation we have as- 
sumed that the incident radiation field deriv es entirely from 
the d ominant ionization star 6 1 Ori C. As IBautista et al.l 
( 1996), we have calculated a dilution factor assuming a T e // 
39000 K, R s tar ~ 9-O Rp) (see j|5^3JI and a distance to the 



Orion Nebula of 414 pc (|Menten et al.ll2007l ). In Table we 
include the Cr + /H + ratio for the nebular and shock compo- 
nent. This is the first estimation of the Cr + abundance in 
the Orion Nebula. 

Several [Fe n] lines have been detected in our spectra. As in 
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Table 3. Ionic abundances and abundance discrepancy fac- 
tors. 



Nebular Component 



Shock Component 



Ionic Abundances from CELs 







t 2 =0 






t 2 >0 






i 2 =0 




t 2 >0 








7.87 






















N+ 


7 


02±0 


04 


7 


07±0 


05 


7 


35±0 


03 


7 


52±0 


04 


o+ 


8 


00±0 


06 


8 


05±0 


09 


8 


29±0 


06 


8 


48±0 


08 


2+ 


8 


35±0 


0.3 


8 


46±0 


01 


8 


08±0 


03 


8 


43±0 


05 


Ne 2 + 


7 


46±0 


11 


7 


58±0 


12 


7 


13±0 


10 


7 


51±0 


11 


S+ 


5 


50±0 


07 


5 


54±0 


08 


6 


03±0 


04 


6 


22±0 


05 


S 2 + 


6 


90±0 


25 


6 


98±0 


25 


6 


89±0 


22 


7 


16±0 


21 


C1+ 


3 


99±0 


09 


4 


04±0 


10 


1 


52±0 


06 


4 


68±0 


06 


Cl 2 + 


5 


13±0 


04 


5 


23±0 


05 


5 


05±0 


05 


5 


38±0 


06 


Ar 2 + 


6 


30±0 


04 


6 


39±0 


04 


6 


26±0 


05 


6 


56±0 


04 


Ar 3 + 


3 


73±0 


11 


3 


85±0 


12 














Ca+ 














3 


86±0 


07 


4 


03±0 


07 


Cr+ 


2 


88±0 


11 


2 


92±0 


11 


3 


75±0 


07 


3 


91±0 


07 


Fe+ 


5 


18±0 


26 


5 


23±0 


27 


5 


82±0 


03 


6 


01±0 


06 


Fe 2 + 


5 


66±0 


13 


5 


72±0 


1.3 


6 


77±0 


09 


6 


96±0 


10 


Fe 3 + 














5 


87±0 


16 


6 


16±0 


20 


Ni+ 


3 


83±0 


10 


3 


88±0 


11 


1 


78±0 


09 


4 


96±0 


09 


Ni 2 + 


4 


42±0 


11 


4 


47±0 


15 


5 


60±0 


09 


5 


77±0 


09 



Ionic Abundances from RLs 



He+ 


10.94±0.01 




10.93±0.01 


C 2 + 


8.32±0.07 




8.25±0.08 


o+ 


8.01±0.12 




8.25±0.16 


o 2 + 


8.46±0.03 




8.44±0.03 


ADFs 


C 2 + 


0.45 






o+ 


0.01±0.17 




-0.04±0.14 


o 2 + 


0.11±0.04 




0.35±0.05 



a In units of 12+log(X+ n /H+) . 

b Av erage value from positions 8b and 11 of Walt er" st al.l 
l|l992h . 



the case of [Cr n] lines , most of them are affected by contin - 
uum fluorescence (see iRodrfguez! Il999l : IVerner et al1l2000h . 
Following the same procedure as for Cr + , we considered a 
159 model atom in order to compute the Fe + abundances 
using the atomic data presented in iBautista fe Pradhanl 
(|l99St ). 

Many [Fe ill] lines have been detected in the two kinematic 
components and their flux is not affected by fluorescence. 
For the calculations of the Fe 2+ /H + ratio, we have imple- 
mented a 3 4-level model atom that uses collision strengths 
ta ken from Zhand J 19961 ) and the transition probabilities 
of iQuinetJ dl996h as w ell as the new transitions found by 
Ijohansson et al.l (|2000r ) . The average value of the Fe 2+ abun- 
dance has been obtained from 31 and 12 individual emission 
lines for the shock and nebular component, respectively. 
One [Fe iv] line has been detected in the shock component at 
6740 A. The Fe 3+ / H + ratio has been derived using a 33-level 
m odel atom where all collisio n strengths are those calculated 
bv lZhang fe Pradhanl dl997h and the transition probabilities 
recommended bv lFroese Fischer fe Rubinl (|2004l ). 
Several [Ni n] lines have been measured in both kinematic 
components but they are strongly affected by continuum flu- 



orescence fsee lLucvl[l995h . As for the Cr + and Fe + ions, we 
have used a 76-level model that includes contin uum fluo- 
rescen ce excitation and the new collisional data of iBautistal 
(2004) in order to compute the Ki + abundances. 
We have measured several [Ni ill] lines in the shock and neb- 
ular components. These lines are not expected to be affected 
by fluorescence. The Ni 2+ /H + ratio has been deri ved using 
a 126 -level model atom and the atomic data of IBautistal 
(|200ll ). 

The final adopted values of the ionic abundances are listed 
in columns (1) and (3) of Table|3]for the nebular and shock 
component, respectively. Columns (2) and (4) correspond 
to the ionic abundances of both components assuming the 
presence of temperature fluctuations (see i]5.5[) . In this table, 
we have also included the C 2 + /H + ratio obtained from UV 
CELs by I Walter et~aH (ll992T ). We have taken the average of 
the values corresponding to their slit positions 8b and 11, 
the nearest positions to HH 202. The uncertainties shown 
in the table are the quadratic sum of the independent con- 
tributions of the error in the density, temperature and line 
fluxes. 

The abundance determinations presented in Table [3] show 
the following behaviour: ionic abundances determined from 
CELs of once ionized species are always higher in the shock 
than in the nebular component; the twice ionized species of 
elements lighter than Ne (included) , show lower abundances 
in the shock than in the nebular component; and the twice 
ionized species of elements heavier than Ne show similar 
abundances in both components except for iron, chromium 
and nickel abundances that show substantially larger abun- 
dances in the shock component, something that can be ex- 
plained if a significant dust destruction occurs in this com- 
ponent (see t|5.6|l . 

Finally, we have compared our abundance determinations 
from UVES data with those obtained in the HH 202-S re- 
gion from the IFU data presented in Paper I. Integrating 
the spaxels of the section of the FOV indicated in i|3] we ob- 
tain 12+log(0 2+ /H+) = 8.18±0.07 from CELs, which is in 
good agreement with the numbers presented in Table [3] for 
both kinematic components. On the other hand, the average 
value of the + abundance from CELs is 8.06±0.14. How- 
ever, considering the large density dependence of this ionic 
abundance, we have recalculated the + /H + ratio adopting 
the physical conditions measured in the shock component 
from UVES, finding a value of 8.26±0.09, which is in bet- 
ter agreement with the UVES determinations for the shock 
component, the brightest one in the [O n] line emission. 

4.4 Ionic abundances from RLs 

We have measured several He I emission lines in the spectra 
of HH 202, both in the nebular and in the shock compo- 
nents. These lines arise mainly from recombination, but they 
can be affected by collisional excitation and self-absorption 
effects. We have used the effecti ve recombination coeffi- 
cients of IStorev fe Hummerl j 19951) for H I and those com- 
puted by iPorter et alj (|2005I)~ with t he interpolation for- 
mulae provided by IPorter et alj l|2007|) for He I. The colli- 
sional contri bution was estimated from Sawev fe Berringtonl 
i| 19931 ) and iKingdon fe Ferlandl i| 19951 ). and the optical 
depth in the triplet line s was derived from the computa- 
tions by iBeniamin et alj (|2002l ). We have determined the 
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He + /H + ratio from a max imum likelihood method (MLM, 
IPeimbert et~aill2000l , l20()3 l 

To self-consistently determine n e (He i), T e (He i), He + /H + 
and the optical depth in the He I 3889 line, T3889, we have 
used the adopted density obtained from the CEL ratios for 
each component as n e (He 1) (see Table [2} and a set of 16 
He 1 lines (at 3614, 3819, 3889, 3965, 4026, 4121, 4388, 4471, 
4713, 4922, 5016, 5048, 5876, 6678, 7065 and 7281 A). We 
have discarded the He I 5048 A line in the nebular compo- 
nent because it is affected by charge transfer in the CCD. 
So, for the nebular component of HH 202, we have a total 
of 16 observational constraints (15 lines plus n e ), and for 
the shock component, we have 17 observational constraints 
(16 lines plus n e ). Finally, we have obtained the best value 
for the 3 unknowns and t 2 by minimizing \ 2 . The final \ 2 
parameters we have obtained are 7.53 for the nebular com- 
ponent and 12.34 for the shock component, which indicate 
very good fits, taking into account the degrees of freedom. 
The final adopted value of the He + /H + ratio for each com- 
ponent is included in Table [3] 

We have detected C 11 lines of multiplets 2, 3, 4, 6 and 17.02. 
The brightest of these lines is C 11 4267 A, which belongs 
to multiplet 6 and can be used to derive a proper C 2+ /H + 
ratio. The rest of the C 11 li nes are affecte d by fluorescence, 
like multiplets 2, 3 and 4 (see lGrandill 19761 ) or are very weak, 
as in the case of the line multiplet 17.02 that has an uncer- 
tainty of 40% in the line flux. 

We have derived the + /H + ratio from RLs for the shock 
and nebular components. The O I lines of multiplet 1 are 
very weak and they are partially blended with bright tel- 
luric emission. In order to obtain the best possible abun- 
dance determination, we have used different lines for each 
component: O I 7775 A for the nebular component and O I 
7772 A for the shock one, these are precisely the lines least 
affected by line-blending. 

The high signal-to-noise of the spectra allowed us to detect 
and measure seven lines of the multiplet 1 of O 11 as we can 
see in Figure [3] These lines are affected by non- local ther- 
mal equilibrium (NLTE) effects (jRuiz et al']l2003l ). therefore 
to obtain a correct 2+ abundance it is necessary to ob- 
serve the eight lines of the multiplet. However, these effects 
are rather small in the Orion Nebula -as well as in the ob- 
served components- due to its relatively large density. Then, 
assuming LTE, the 2+ abundance from RLs has been cal- 
culated considering the abundances obtained from the flux 
of each line of multiplet 1 and the ab undance from the est i- 
mated total flux of the multiplet (see lEsteban et afl ll998). 
The abundance determinations in Table [3] show that the 
He+/H+, C 2+ /H+ and 2+ /H+ ratios derived from RLs are 
always very similar in both shock and nebular component. In 
the case of + abundances, the nominal values determined 
for both components seem to be somewhat different (about 
0.24 dex), but they are marginally in agreement considering 
the large uncertainties of this ion abundance. 
As in the previous subsection, we have compared our abun- 
dance determinations from RLs with those obtained for 
HH 202-S in Paper I. From the data of Paper I, we ob- 
tain 12+log(0 2+ /H+) = 8.39L0.13 and 12+log(C 2+ /H+) = 
8. 29L0.11, values which are in good agreement with those 
obtained in this paper. 



4.5 Abundance discrepancy factors 

We have calculated ionic abundances from two kinds of lines 
-RLs and CELs- for three ions: C 2+ , 0+ and 2+ . We 
present their values for the two components in Table [3] We 
have computed the ADF for these ions using the following 
definition: 



ADF(X-)=log^j RL -log(g: /c]L 



(3) 



In the case of the ADF(C 2+ ) it can only be estimated for 
the nebular component and from the comparison of our 
determination from RLs and those from CELs for nearby 
zones taken from the literature. The value of the ADF(C 2+ ) 
amounts to 0.45 dex. 

On the one hand, the ADF(0 + ) can also be estimated in 
our spectrum and shows values very close to zero. However, 
these ADF values are rather uncertain. On the other hand, 
as we can see in Table [3] the 2+ abundance from RLs is 
the same for both components while that from CELs is lower 
in the shock component, probably because the recombina- 
tion rate increases in the shock one. This fact produces an 
ADF(0 2+ ) about 0.2 dex higher in the shock component 
than in the nebular one. This striking result will be dis- 
cussed in section £| 5 . 5 1 

The values of the ADFs of C 2+ and 2+ for the nebular 
component ar e in go od agreement with those obtained by 
lEsteban et ail (|2004l ) for a zone closer to the Trapezium 
cluster than HH 202. I n the case of the ADF (Q+), both de- 
terminations disagree, lEsteban et ail (|2004 ) report a much 
larger value (0.39 dex). 

4.6 Total abundances 

In order to derive the total gaseous abundances of the dif- 
ferent elements present in our spectrum, we have to correct 
for the unseen ionization stages by using a set of ionization 
correction factors (ICFs). The adopted ICF values are pre- 
sented in Table [3] and the total abundances in Table [S] As 
in the case of the ionic abundances from CELs, these tables 
include values under the assumption of t 2 = -columns (1) 
and (3)- and under the presence of temperature fluctuations 
(see Sj53} -columns (2) and (4)-. 

The total helium abundance has been corrected for the pres- 
ence of neutra l heliu m using the expression proposed by 
IPeimbert etafl i|l992l ) based on the similarity of the ion- 
ization potentials of He (24.6 eV) and S+ (23.3 eV). 



He 

IT 



He + T ™ /TT on He + 
x -=z— = ICF(He ) x 



(4) 



For C we have ado pted the ICF(C + ) der ived from photoion- 
ization models of IGarnett et al.l (jj.999) for the shock and 
nebular component. In order to derive the total abundance 
of nitrogen we have used the usual ICF: 



N 
H 



Q+ +Q 2 
0+ 



x = ICF(N 2+ ) 



N+ 

x w . (5) 



This expression gives very different values of the ICF(N 2+ ) 
for both components due to their rather different ionization 
degree. 

The total abundance of oxygen is calculated as the sum of 
+ and 2+ abundances. The absence of He 11 lines in the 
spectra and the similarity between the ionization potentials 
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Table 4. Adopted ICF values. 





Nebular Component 


Shock Component 


Elements Unseen Ion 


t 2 =0 t 2 >0 


t 2 =0 t 2 >0 



He 


He" 


1.04±0.02 


1.03±0.02 


1.12± 0.06 


1.08±0.04 


C a 


C+ 


1.31±0.46 


1.50±0.47 


N 


N 2 + 


3.21±0.54 


3.82±0.83 


1.62± 0.27 


2.29±0.46 


Ne 


Ne+ 


1.45±0.15 


1.35±0.19 


2.61± 0.31 


1.78±0.26 


S 


g3+ 


1.01±0.01 


l.OliO.Ol 


1.09± 0.03 


1.03±0.01 


Ar b 


Ar+ 


1.33 








Ar c 


Ar+ 


1.20±0.36 


1.16±0.36 


2.00± 0.51 


1.71±0.46 


Fc d 


Fe 3 + 


2.71±0.46 


3.16±0.69 


1.52± 0.25 


2.02±0.41 


Ni 


Ni 3 + 


3.21±0.54 


3.82±0.83 


1.62± 0.27 


2.29±0.46 



a From photoionization models bv lGarnett et aL | jl999T) . 
b Mean of Orion Nebula models. 

c From correlations obtained by [Martin-Hernandez et al.l d2002T ). 
d From photoionization models by Rodriguez & Rubin (2005). 



of He + and 2+ , implies the absence of O s+ . In Table [S] we 
present the O abundances from RLs and CELs. 
The only measurable CELs of Ne in the optical range 
are those of Ne 2+ but the fraction of Ne + can be impor- 
tant in the nebula. We hav e adopted the usual expression 
l|Peimbert fe Costerol Il969l ) to obtain the total Ne abun- 
dance: 

+ + 2+ Ne 2+ j_ Ne 2+ 

ICF(Ne+)x^^. (6) 



Ne 

IT 



o 2 + 



H^ 



We have measured CELs of two ionization stages of S: S + 
and S 2+ . Then, we have used an IC F to take into account 
the presence of S 3+ (|Stasinskall 19781 ) which is based on pho- 
toionization models of H n regions. 



H 



H 



l o+ + o 2 +- 



-i/3 S + + S 
x 



2 + 



rL 



= ICF(S 3+ ) x ^-±S 



Following lEsteban et alj (|l998l ) , we expect that the amount 
of Cl 3+ is negligible in the Orion Nebula. Therefore, the to- 
tal abundance of chlorine is simply the sum of C1+ and Cl 2+ 
abundances. 

For argon, we have determinations of Ar 2+ and Ar 3+ 
but some contribution of Ar + is expected. In Table [3] 
we present the values obtained from two ICF schemes: 
one obtained from correlations between N 2+ /N + vs. 
Ar 2+ /Ar + from ISO observatio ns of compact H II regions by 
Martfn-Hernandez et al.l (|2002t ) and another one -following 
Osterbrock et alj 1119921)- derived as t he mean of Orion Neb - 
ula models bv lRubin et all (|l99ll ) and lBaldwin et alj i|l99ll ). 
We have measured lines of three ionization stages of iron in 
the shock component -Fe + , Fe 2+ and Fe 3+ - and two stages 
of ionization in the nebular component -Fe + and Fe 2+ -. For 
the shock component, we can derive the total Fe abundance 
from the sum of the three ionization stages. For the nebular 
component -and also for the shock one in order to compare- 
we ha ve used an ICF scheme base d on photoionization mod- 
els of iRodrfguez fe Rubinl (|2005l ^ to obtain the total Fe/H 
ratio using only the Fe 2+ abundances, which is given by: 



Fe / + 



cos Fe 2+ O 
x - , x 



0+ 



H 



(8) 



Finally, there is no ICF available in the literature to cor- 
rect for the presence of Ni 3+ in order to calculate the total 
Ni abundance. Nevertheless, we have applied a first order 
ICF scheme based on the similarity between the ionization 
potentials of Ni 3+ (35.17 eV) and 2+ (35.12 eV): 



Ni 3 



Ni 



Therefore, 



Q 2 + 
O 

Ni 4 



Ni _ _0_ /Ni+ Ni 2+ \ 
~R ~ 0+ X \B+ + ~W~)' 



(9) 



(10) 



In general, the total abundances shown in columns (1) and 
(3) of Table [5] are quite similar for the shock and nebular 
component within the errors, except for the nickel and iron 
abundances, which are much larger in the shock component 
-see section t]5.6l for a possible explanation. The set of abun- 
dances for the nebular com ponent are in very go od agree- 
ment with previous results of lEsteban et al.1 (|2004l ') . We have 
also compared o ur Ni abundance values with the previous 
determination of lOsterbrock et al.1 (|l992t l finding that our 
Ni/H ratio for the nebular component is an order of magni- 
tud lower. This difference is due to the larg e uncertainties in 
the atomic data used by those authors (seefBautista 200lJ ). 



5 DISCUSSION 

5.1 Differences between the c(H/3) coefficient 
determined with different lines 

A puzzling feature of our UVES data is that the c(H/3) values 
determined with different lines ratios appear to be inconsis- 
tent with each other, even with observational errors taken 
into account. Possible explanations are either a bias in the 
extinction curve or an extra mechanism, in addition to ex- 
tinction, altering the individual line intensities from case B 
predictions. 

Figure [4] shows the c(H/3) values measured for individual 
line ratios in the Balmer and Paschen series. The bizarre 
pattern followed by the curve, and particularly the steep 
slope it reaches in the proximity of the Balmer and Paschen 
limits, strongly suggest that the solution cannot be a bias 
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Table 5. Total abundances" 



Nebular Component 



Shock Component 





t 2 =0 


i 2 >0 


t 2 =0 


t 2 >0 


Ho 
XI c 


10.95±0.01 


10.95±0.02 


10.98±0.03 


10.98±0.02 


c b 


8.07 








c c 


8.43±0.17 


8.43±0.16 


N 


7.53±0.08 


7.62±0.11 


7.56±0.08 


7.81±0.10 


o 


8.51±0.03 


8.60±0.04 


8.50±0.04 


8.76±0.05 


o c 


8.59±0.05 


8.65±0.05 


No 


7.62±0.12 


7.72±0.13 


7.54±0.11 


7.83±0.12 


S 


6.92±0.24 


7.00±0.24 


6.98±0.19 


7.23±0.19 


CI 


5.16±0.04 


5.26±0.05 


5.16±0.04 


5.46±0.05 


Ar d 


6.42±0.04 


6.52±0.04 






Ar e 


6.38±0.19 


6.46±0.14 


6.56±0.21 


6.79±0.12 


Fef 






6.86±0.07 


7.06±0.08 


p e s 


6.10±0.15 


6.19±0.16 


6.95±0.12 


7.19±0.13 


Ni 


5.03±0.14 


5.12±0.15 


5.87±0.11 


6.11±0.12 



a In units of 12+log(X+ n /H+) . 

b Av erage value from positions 8b and 11 of Walt er" st al.l 
jl992h . 
c Value derived from RLs. 

d Adopting the ICF from the mean of Orion Nebula models. 
e Adopting the ICF from ISO observations 
jMartm-Hernandez et al]|2002h . 
f From Fe++Fe 2 ++Fe 3 +. 
9 Assuming the ICF of equation (6). 

in the extinction law. 

We can test the second hypothesis by considering how the 
line intensity ratios relative to case B ratios depend on 
the principal quantum number n of the level where each 
line originates. Such dependence, plotted in Figure [5] for 
the two kinematic components, shows a definite trend with 
n for both series. This strongly supports the second of 
our hypotheses, namely that an extra mechanism is act- 
ing to deviate level populations away from case B predic- 
tions. Indeed, detailed photoionization modeling indicates 
that this behaviour is the result of two independent but 
concomitant mechanisms: ^-changing collisions with |A£| > 
1 and pumping of Balmer and Paschen lines by absorp- 
tion of stellar continuum photons at the Lyman wave- 
lengths. Botii_jnechanisms are neglected in case B calcu- 
lations l|Storev fc Hummerlll995h but included in our mod- 
els, which take advantage of a new model hydrogen atom 
with fully resolved levels (Cloudy, version C08.00; Porter, 
Ferland, van Hoof, & Williams, in preparation; see also Ap- 
pendix A in Luridiana et al. 2009); both alter the n, I popu- 
lations, resulting in enhanced intensities of the high n lines. 
As for the first mechanism, at low n it has a negligible ef- 
fect if compared to other depopulation mechanisms, such 
as energy-changing collisions and horizontal collisions with 
A£| = 1; at high n, it becomes increasingly important. Case 
B calculations neglect this mechanism by construction, so a 
discrepancy is doomed to appear whenever high-n lines are 
compared to case B results. 

The effectiveness of the second mechanism strongly de- 
pends on the availability of the exciting photons, i.e. on 
the stellar flux at the Lyman wavelengths; the results of 
iLuridiana et all l|2009h and further preliminary calculations 
(Luridiana et al. in preparation) suggest that its impact on 
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Figure 4. c(H/3) values obtained from different line ratios vs. 
the inverse of the wavelength in /im for the nebular (up) and 
shock component (down). The filled and non-filled squares cor- 
respond to the c(H/3) values derived from Balmer and Paschen 
lines, respectively. The horizontal line is the weighted average 
value adopted and the grey band its error 
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Figure 5. Dereddened fluxes of Balmer (filled squares) and 
Paschen (non-filled squares) lines to their theoretical flux under 
the Case B prediction ratio vs. the principal quantum number n 
for the nebular (up) and shock component (down). 



line intensities might increase with n. 

A full account of b o th pro cesses in H II regions can be found 
in ILuridiana et al] (|2009l ) and Luridiana et al. (in prepara- 
tion). 



5.2 Comparison of line ratios in the shock and 
nebular components and the ionization 
structure 

To maximize the shock-to-nebula ratio, the echelle spectra 
were extracted over the area where the shock component 
is brighter and the velocity separation with respect to the 
nebular background gas is maximum (see Figure [3). 
In Figure |6j we present the weighted average shock-to- 
nebular ratio for different ionic species, I(A) s h//(A) no b - 
which was defined in equation (1)- with respect to the ion- 
ization potential (IP) needed to create the associated origi- 
nating ion. In general, as we can see in this figure, the line 
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the bulk of the 2+ emission is located at the south of the 
extracted area, O at the north and + is located approxi- 
mately at the centre of the extracted area. Another interest- 
ing feature that can be seen in Figure[2]is that the emissions 
from RLs and CELs of the same ion seem to show the same 
spatio-kinematic profiles. These profiles provide further in- 
dication that HH 202-S is an ionization-bounded shock. 
The peak emission of the spatio-kinematic profiles of ions 
with IP lower or similar to the one of + : Cr + , Ni + , S + , 
Fe + or H+ is located to the north of HH 202-S, as in the 
case of + . Ions as He + , C 2+ or Fe 2+ show their peak emis- 
sion about the centre of the aperture, while the ions with the 
highest IP, as Ne 2+ or Ar 3+ , show spatio-kinematic profiles 
similar to that of 2+ . 



Figure 6. Weighted average of shock-to-nebular line ratios vs. 
the ionization potencial (IP) needed to create the associated orig- 
inating ion. Non-filled circles correspond to ratios measured from 
RLs (only for O+ and 0+ + ). Smaller error bars are of about 3%. 

ratios of the shock component relative to those of the ambi- 
ent gas are between 1 and 2 for most ionic species with IP 
above 10 eV, and close to one for the most ionized species as 
2+ and Ne 2+ . Since the illumination of the shock should 
be approximately the same as that of the nebula at this par- 
ticular zone of the slit, shock-to-nebular ratios of the order 
of one imply that the shock should be ioniz ation-bounded. 
This i s exactly the opposite situation that iBlagrave et al.l 
(2006) find in HH 529, where the shock-to-nebular ratios 
are clearly lower than 1 indicating that the shock associated 
to HH 529 is matter-bounded. 

As we have commented above and can be seen in Figure [5] 
the shock-to-nebular ratio varies from 1 to 2 for species with 
an IP higher than 10 eV except in two cases: a) Fe 2+ , whose 
lines show a substantial enhancement in the shock compo- 
nent, probably due to dust destruction (see N5.6p . and b) 
Cl + , but this can be accidental because this ionic abundance 
is derived from a single rather faint line with a very large 
uncertainty. Ionic species with an IP lower than 10 eV show 
a shock-to-nebular ratio always higher than 2 in Figure [6] 
In particular, Fe + , Ni + and Cr + show ratios larger than 2 
and these ions may be also affected by an increase of the 
gas-phase abundance due to dust destruction. 
Neutral species like O or N° are associated with the pres- 
ence of an ionization front. In the case of O , the shock-to- 
nebular ratio has been calculated from the [O i] 6300 and 
6363 A lines that are contaminated by telluric emissions but 
only in the shock component. We extracted the telluric emis- 
sions from the zone free of shock emission along the slit and 
subtracted this feature to the shock component. The shock- 
to-nebular ratio of [O i] lines is the largest one for those ions 
which are not heavily affected by possible dust destruction. 
This is a further indication of the presence of an ionization 
front in HH 202-S. In the case of [N i] lines, those belonging 
to the nebular component are also contaminated by telluric 
emission but, in this case, it was impossible to deblend prop- 
erly these lines. 

In Figure [2] we can see the spatio-kinematic profiles of lines 
of different ionization stages of oxygen: O , + and 2+ , 
including RLs and CELs. In the shock component, we can 
clearly distinguish a stratification in the location of the ions: 



5.3 Width of the ionized slab and the physical 
separation between (^Ori C and HH 202-S 

An interesting result of this paper is the claim that HH 202- 
S contains an ionization front as we have shown in >14.1l and 
£15.21 Due to this fact, we can estimate the width of the 
ionized slab of HH 202-S and its physical separation with 
respect to the main ionization source of the Trapezium clus- 
ter, 6 1 Ori C. 

On the one hand, from the maximum emission of the shock 
component in the spatio-kinematic profiles of 2+ and O 
shown in Figure [2] we can measure an angular distance on 
the plane of the sky of about 3'. / 9±0'. , 5. Using t he dis tance 
to the Orion Nebula obtained by iMenten et al.1 ((2007), d = 
414±7 pc, and the inclination angle of HH 202-S with re- 
spect to the plane of the sky calculated bv lO'Dell fc Hennevl 
<|200Sl ). 9 = 48°, we estimate (11.7±1.5) x 10" 3 pc for the 
width of the ionized slab. 

On the other hand, to trap the ionization front in HH 202-S, 
the incident Lyman continuum flux must be balanced by the 
recombinations in the ionized slab, i.e., 

FLy = 1^ = n sh a b(H°,T)L, (11) 

where D is the physical separation between HH 202-S and 
9 1 Ori C, Q(H°) is the ionizing photon rate, n s h the density 
in the shock component, as(-ff°, T) is the case B recombina- 
tion coefficient for H and L the width of the slab. In order to 
estimate Q(H°), we have used a spectral energy distribution 
of fastwin d code with the stellar p arameters for 9 1 Ovi C 
obtained bv lSimon-Dfaz et al.1 l|2006r ) -T eff = 39000 ± 1000 
K and log g = 4.1 ± 0.1 dex- and the dis tance to the Orion 
Nebula calculated bv lMenten et al.1 (120071 ). Then, the output 
parameters have been: the stellar radius R = (9.0 ± 1.3)-R0, 
the stellar luminosity log L — 38.80 ± 0.14 dex and the ion- 
izing photon rate Q(H°) = (6.30 ± 2.00) x 10 48 s _1 . Taking 
a value for the recombination coefficient, as = 2.59 x 10 -13 
cm " 3 s" 1 at 10 4 K, we have finally calculated a physical 
separation of D — 0.14 ± 0.05 pc. This result suggests that 
HH 202 is quite embedded within the body of the Orion Neb- 
ula and, therefore, discarts the origin of the ionized front as 
result of the interaction of the gas flow with the veil (see 
i]4.ip , whic h is between 1 and 3 pc in front of the Trapezium 
cluster (see lAbel et al.ll2004l ). 
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Figure 7. Heliocentric velocities of the nebular (up) and shock 
(down) component emission lines us. the ionization potencial (IP) 
needed to create the associated originating ion. The velocity of the 
photon dominated region (dashed horizontal line) is also shown. 
Non-filled triangles correspond to ions observed in only one kine- 
matic component. Non-filled circles correspond to heliocentric ve- 
locities measured from RLs (only for + and +_l "). 



5.4 Radial velocity analysis 

In Figure [7] we show the average heliocentric velocity of 
the lines that belong to a given ionic species as a function 
of the IP needed to create the associated originating ion. 
We have used two separate graphs to distinguish between 
the behaviour of the nebular and the shock component. 
The nebular component shows the typical velocity gradient 
that has been observed in other positions of the Orion Neb- 
ula and other H 1 1 region s (e.g. iBautista fc Pradhar]|l998l ; 
lEsteban fc Peim bcrt 1999]), with a velocity difference of the 
order of —15 km s _1 between the neutral and the most ion- 
ized species. This gradient is likely produced by the presence 
of flows of ionized gas originating from the ionization front 
inside the nebula. In contrast, the ions in the shock compo- 
nent show a rather similar radial velocity independently of 
their IP. This indicates that the bulk of the ionized gas at 
the HH 202-S is moving at approximately the same velocity 
with respect to the rest of the nebula. The H I lines (Balmer 
and Paschen series) of the shock component are shifted by 
— 50.7±1.0 km s -1 relative to the H I lines of the nebular 
component and — 64.8±1.0 km s~ relative to the velocity 
of the photon dominated region (PDR, +28 km s- MCoudis! 
1 19821 ). The zone covered by our UVES slit ( see F igure ED) 
coincides with position 117-256 of lDoi et af] l|2004 ). These 
authors detect two radial velocity components belonging to 
the shock gas in this zone, a fast component ( — 57 km s _1 ) 
and a slower brighter one (—31 km s~ ). Our value of the ra- 
dial velocity for the shock component of HH 202 (— 36.8±1.0 
km s _1 ) is so mewhat more neg ative than the slower velocity 
component of iDoi et alj (|2004| ). probably our shock compo- 
nent corresponds to the unresolved blend of the two velocity 
systems detected by those authors. 



5.5 The abundance discrepancy and temperature 
fluctuations 

Assuming the validity of the temperature fluctuations 
paradigm and that this phenomenon produces the abun- 



dance discrepancy, we can estimate the values of the t pa- 
rameter from the ADFs obtained for each component and 
ion (see Table [3]). In Table HJ] we include the t 2 values that 
produce the agreement between the abundance determina- 
tions obtained from CELs and RLs of 2+ and C 2+ . These 
calculati ons have been made follow ing the formalism out- 
lined bv lPeimbert fc Costerd (1 19691 ). Adopting the t 2 value 
obtained for 2+ zone, t 2 (0 2+ ), we have calculated the ionic 
abundances, the ICFs and the total abundances under the 
presence of temperature fluctuations and they are presented 
in Tables [3] [4] and [5] respectively. 

As we can see in Table [S] the high t 2 (0 2+ ) value used to de- 
rive the abundances for the shock component implies that 
the total abundances obtained from CELs considering t 2 >0 
are higher in the shock component than in the nebular one 
for all the elements. On the other hand, the abundances of 
C and O obtained from RLs are very similar in both com- 
ponents. Moreover, the O/H ratio computed from CELs for 
t 2 >0 is higher than that obtained from RLs. However, for 
the nebular component, we find that the total abundance of 
oxygen determined from CELs considering t 2 >0 agrees with 
the oxygen abundance from RLs. 

These results can suggest that, perhaps, the t 2 value we have 
found for the shock component is too high and, therefore, 
that the t 2 paradigm is not applicable in this case. However, 
we have to consider that we have used a t 2 value representa- 
tive of the high ionization zone for all the ions. The increase 
of the + /0 2+ ratio in the shock component with respect 
to the nebular one, due to an enhanced recombination rate, 
makes the + zone more extended in this component and a 
lower t 2 (0 + ) would lead to a better agreement between the 
O abundances of both components. We have also computed 
the t 2 parameter for the He + zone using a MLM (see H4.4[) . 
The determination of the i 2 (He + ) weighs the + and the 
2+ zones depending on their extension. It is remarkable 
that the t 2 values obtained from different methods, such as 
the He + lines and the ADF(0 2+ ), which assume the pres- 
ence of temperature fluctuations in the observed volume, 
produce almost identical values in both components, though 
t 2 (He + ) has large uncertainties. These ones can reconcile the 
total abundances of both components for t 2 >0 considering 
that the MLM depends strongly on the He I 3889 A line 
flux and that the shock components of this line and H8 are 
severely blended. 

In Table EH we have also included the values of t 2 obtained 
bv lBlaerave et al.l l|2006l ) from the ADF(0 2+ ) that they esti- 
mate for the nebul ar and shock compon ents of HH 529 and 
those obtained by lEsteban et all (l2004h from the ADF of 
2+ and C 2+ . The data of lEsteban et"al] (|2004l ) correspond 
to a zone free of high velocity flows and can be considered 
as representative of the nebular component but closer to the 
Trapezium. The values of the t 2 obtained for the nebular 
compone nt of HH 202-S are qu ite consistent with those ob- 
tained by lEsteban et al" | i|2004l ). T he t 2 (Q 2+ ) of th e nebu lar 
component of HH 529 obtained by iBlagrave et al.l (2006) is 
lower than the two other determinations, but still consistent 
with our values within the uncertainties. As has been com- 
mented above, the t 2 values found in the shock component 
of HH 202 are much higher than tho se of the nebular com - 
ponent and the value determined by IBlagrave et all (|2006D 
for the shock component of HH 529. 

The effect of temperature fluctuations in the spectra of ion- 
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Table 6. Estimates of the t 2 parameter. 





HH 202 


HH 529 a 






Nebular 


Shock 


Nebular 


Shock 




Ion 


Component 


Component 


Component 


Component 


Esteban et al. ('2004') 


o 2 + 


0.016±0.006 


0.050±0.007 


0.009±0.004 


O.OlOiO.010 


0.022±0.002 


C 2+ 


0.040 








0.039±0.011 


He+ 


0.014±0.013 


0.049±0.011 









a Blagrave etafl (2006). 



ized nebulae and their existen ce has been a controversial 
problem from the first work of Peimbert] J 19671 ). Peimberj 
i| 19951 ). lEstebanl (|2002l ). iPeimbert fe PeimbertJ (|2006l ) have 
reviewed the possible mechanisms that could produce such 
temperature fluctuations. Among the possible sources of 
temperature fluctuations, we can list the deposition of me- 
chanics ene rgy by shocks and th e gas compresion due to 
turbulence. IPeimbert et al.l |l991) have studied the effect 
of shock wave s in H II regions using shock models by 
lHartigan et al.l (1 19871 ). They found that high t 2 values can 
be explained by the presence of shocks with velocities larger 
than 100 km s _1 . In the case of the flow of 89 km s _1 that 
produces HH 202, this effect should produce very small t 2 
values. Therefore, another process or other formalisms of 
the problem should be taken into account in this case. Fur- 
ther duscussion of this problem is beyond the scope of this 
article. 



5.6 Dust destruction 



As indicated in H5.2I the emission lines of refractory ele- 
ments as Fe, Ni or Cr are much brighter in the gas flow 
than in the ambient gas. Moreover, the shock component 
shows relatively bright [Ca n] lines which are not detected 
in the spectrum of the ambient gas. On the other hand, in 
i]4.1l we have stated that HH 202-S does not contain a sub- 
stantial contribution of shock excitation and, therefore, we 
need other mechanisms to explain those abnormally high line 
fluxes. It is well known that Fe, Ni, Cr and Ca are expected 
to be largely depleted in neutral and molecular interstellar 
clouds as well as in H II regions. However, theoretical studies 
have shown that fast shocks -as those typical in HH objects 
and supernova remnants- should efficiently destroy grains 
by thermal and nonthermal sputtering in th e gas behind 
the shock front and by grain-grain collisions (IMcKee et all 
1 1984 I Jones et al J 1 1994 iMouri fe Taniguchl |2000| ) . Several 
works have shown that some nonphotoionized HH flows show 
a decrease in the amount of Fe depletion as determined 
from the analysis of [Fe h] lines ( Beck-Winchatz et al.ll 19961 ; 



iBohm fe Mattll200ll ; iNisini et al 



20051 ). On the other hand, 



HH 399 is the first precedent of a fully ionized HH object 
where an overabundance in Fe w as detected and this one 
was related with dust destruction (Rodriguez 2002). 
In Table[7]we compare the values of the C, O, Fe, Ni and Cr + 
abun dances and the Fe/Ni ratios of the Sun (|Grevesse et alj 
120071 ) and the nebular and shock components for t 2 =0 and 
t 2 >0. In the case of Fe and Ni, we can see that the difference 
between the abundances of the shock and nebular compo- 
nents are of the order of 0.85 dex for t 2 —0 and 1.00 dex for 
t 2 >0. This result indicates that the gas-phase abundance of 



Fe and Ni increases by a factor between 7 and 10 after the 
passage of the shock wave. The fact that the Fe/Ni ratio 
is the same in both components -because the increase of 
the gas-phase abundance of both elements is the same- and 
consistent with the solar Fe/Ni ratio within the errors, sug- 
gests that the abundance pattern we see in HH 202-S is the 
likely product of dust destruction. In fact, observations of 
Galactic interstellar clouds indicate tha t Fe, Ni -as well as 
Cr- have the same dust-phase fraction (|Savage fe Sembachl 
1996b; Jones 2000). It is also remarkable that, although it 
could be modulated by ionization effects, the behaviour of 
the Cr + abundance is also consistent with the dust destruc- 
tion scenario. The increase of the Cr + /H + ratio between the 
nebular and shock components is identical to that of Fe and 
Ni. 

Considering the solar F e/H ratio as the reference 
(12+log(Fe/H) = 7.45+0.05. iGrevesse et aLll2007f l -which is 
almost identical to the Fe/H ratio determined for the B-type 
stars of the Orion as sociation (12+log(Fe/H) = 7.44±0.04, 
iPrzvbilla et al.1 [20081 )- we estimate that the Fe dust-phase 
abundance decreases by about 30% for t 2 =0 and about 53% 
for t 2 >0 after the passage of the shock wave in HH 202- 
S. This result i s in good agreeme nt with the predictions 
of the models of I Jones et al l (|1994| ) and[M ouri fe Taniguchil 
(2000). In p articular, for the velocity determined f or HH 202 
(89 km s- MO'Dell fe Hennevllioolh . lJones etafl (|l994l ) ob- 
tain a level of destruction of iron dust particles of the order 
of 40%. On the other hand, the Fe gas-phase abundance we 
measure in HH 202-S follows closely the empirical correlation 
obtained by IBohm fe Matt] (|200lf ) from [Fe 11] and [Ca 11] 
emission line fluxes of several nonphotoionized HH flows. 
Further evidence that the dust is not completely destroyed 
in HH 202 is the detection of 11.7 /im emissio n coincident 
with the bright ionized gas around HH 202-S |Smith et al.l 
2005). As we can see in Table [7J the depletion factors of 
Fe and Ni in the nebular component are similar to those 
found in warm neutral interstellar environments and those 
at HH 202-S of t he order of the de pletions observed in the 
Galactic halo (see lWeltv et al.lll999l . and references therein). 
In the cases of C and O, the effect of dust destruction in their 
gas-phase abundance is more difficult to estimate. Firstly, 
these elements are far less depleted in dust grains than Fe or 
Ni in neutral interstellar clouds and in H 11 regions. Secondly, 
there is still a controversy a bout the correct s olar abundance 
of th ese two elements (see iHolwegerl 1200 ll ; iGrevesse et al.l 
2007). Finally, chemical evolution models predict some in- 
crease in the C/H and O/H ratios in th e 4.6 Gyr since th e 
formation of the Sun (0.28 and 0.13 dex: ICarigi et al1l2005h . 
All these problems make impossible to estimate confident 
values of the depletion factors for C and O. 
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Table 7. Comparison of abundances and depletion factors 







Nebular Component 


Shock Component 






t — U i J> u 


l — U i J> u 


c 


8.39±0.05 


8.43±0.17 


8.43±0.16 





8.66±0.05 


8.59±0.05 


8.65±0.05 


Fe 


7.45±0.05 


6.10±0.15 6.19±0.16 


6.95±0.12 7.19±0.13 


Ni 


6.23±0.04 


5.03±0.14 5.12±0.15 


5.87±0.11 6.11±0.12 


Fe/Ni 


1.22±0.06 


1.07±0.23 1.07±0.25 


1.08±0.17 1.08±0.19 


Cr+ 




2.88±0.11 2.92±0.11 


3.75±0.07 3.91±0.07 




Shock— Nebular 


Nebular— Sun 


Shock— Sun 




t 2 =0 i 2 >0 


t 2 =0 i 2 > 


i 2 =0 i 2 >0 


C 


0.00±0.27 


0.04±0.19 


0.04±0.18 





0.06±0.07 


-0.07±0.07 


-0.01±0.07 


Fe 


0.85±0.13 1.00±0.14 


-1.35±0.05 -1.26±0.05 


-0.50±0.06 -0.26±0.09 


Ni 


0.84±0.12 0.99±0.13 


-1.20±0.04 -1.11±0.04 


-0.36±0.06 -0.12±0.10 


Fe/Ni 








Cr+ 


0.87±0.07 0.99±0.07 







a iGrevesse et alj j200Tt) . 



Table 8. Input parameters for photoionization models. 



Parameter 


Nebular Model 


Shock Model A 


Shock Model B 


Shock Model C 


log(rie) (cm -3 ) 


3.46 


4.24 


4.24 


4.24 


log(C7) 


-1.95 


-2.53 


-2.53 


-2.53 


Dust 


Orion 


Orion 


Orion 


0.5 X Orion 


(Fe/H) 9a3 


2.1 x icr 6 


2.1 X 10" 6 


1.5 X 10" 5 


1.5 x 10" 5 



Considering the data gathered in Table [7] the C abundance 
is virtually the same in the nebular and shock components, 
although any possible small difference may be washed out 
due to the intrinsic relatively large error of the abundance 
determination of this element. In any case, only a slight in- 
crease of the C gas-phase abundance would be expected after 
the passage of the shock wav e for a shock velocity of about 
89 km s" 1 (jjones et al.lll994l ). 

The uncertainties in the determination of the O abundances 
are much lower than in the case of C and its determina- 
tion does not depend on the selection of an appropriate ICF 
scheme. In Table [7J we see an increase of 0.06 dex in the 
O abundance in the shock component with respect to the 
nebular one. These results suggest a possible moderate de- 
crease of the depletion for this element in the shocked ma- 
terial. Very recent detailed determinations of the O abun- 
dance of B-type stars in the Orion association (Simon-Diaz, 
in preparation) indicate that the mean O/H ratio in this 
zone is 8.76±0.04. In principle, one would expect that this 
is a better reference for estimating the dust depletion in the 
Orion Nebula than the solar one because it corresponds to 
the actual O abundance at the same location. If we take 
the B-type stars determination as a reference for the total 
O/H ratio, the amount of O depletion in the ambient gas 
would be -0.17±0.06 dex and -0.11±0.06 dex in the gas 
flow. Therefore, a 30% of the O tied up into dust grains 
would be destroyed after the passage of the shock front. 
There are two other alternative methods to estimate the 
ox ygen depletion factor. The first one can be drawn follow- 
ing [EstebanelaLj (|l998l). This one is based on the fact that 
Mg, Si and Fe form molecules with O which are trapped 



in dust grains. We have used Mg, Si and Fe depletions in 
order to obtain the fraction of O trapped in dust grains. 
These depletions are estimated con sidering the Or i on ga s 
abundances of Mg and Si given by lEsteban et all |l998), 
our Fe abundance for t 2 >0, the O abundances derived 
from RLs, the stellar abundances of the Orion association 
of Si from Simon-Dia z (7.14±0.10) and Mg and Fe from 
IPrzvbiila" et all (|2008l ). By assuming that O is trapped in 
olivine (Mg, Fe^SiOi, pyroxene (Mg, Fe)SiQ3 and several 
oxides like MgO, Fe203 and Fe304 fsee lSavage fe Sembachl 
Il996al . and references therein), we have estimated an O 
depletion factor of — 0.10±0.04 dex for the nebular com- 
ponent. If we consider the older ste llar abundances of Si 
and Fe from Cunha fe Lambert] (|l994l ') and the Mg/Fe ratio 
from lGrevesse et al. (2007), this value becomes — 0.08±0.05, 
which is almost identical to that computed from the new 
stellar Mg, Si and Fe abundances. 

The last method to obtain O depletion factors assumes that 
oxygen and iron are destroyed in the same fraction. The 
fraction of iron dust particules is calculated using the Fe 
abundances for t 2 >0 presented in Table [7] and Fe abundance 
of B-type stars of the Orion association. From this assump- 
tion, and taking into account the O abundance from RLs, 
we have derived a depletion factor of — 0.11 J^ 11 d ex m the 
nebular component. Finally, we have calculated the average 
of the values obtained from the three methods, finding that 
the O depletion factor of the ambient gas is — 0.12±0.03. In 
all cases, the depletion becomes larger by about 0.08 dex 
if we adopt abundances for t 2 =0. Aditionally, the depletion 
factor in the shock component is smaller than in the ambient 
gas, probably due to gas destruction. 
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Table 9. Results of photoionization models. 



Constraint 


Nebular Obs. 


Nebular Mod. 


Shock Obs. 


Shock Model A 


Shock Model B 


Shock Model C 


n e (cm" 3 ) 


2890±550 


2960 


17430±2360 


17540 


17670 


17650 


T e ([N II]) (K) 


9610±390 


9340 


9240±300 


9510 


9220 


9210 


T e ([0 III]) (K) 


8180±200 


8230 


8770±240 


8810 


8670 


8640 


log(0+/0 2 +) 


-0.35±0.07 


-0.37 


0.20±0.07 


0.21 


0.23 


0.22 


J([Fe m] 4658)//(H/3) 


0.009±0.001 


0.010 


0.110±0.006 


0.017 


0.111 


0.107 



5.7 Photoionization models for HH 202-S 

To test the role of dust destruction both in the temperature 
structure observed in the shock component and in the iron 
abundance in the gas, we have run some simple photoioniza- 
tion models for the nebular and the sh ock components. Th e 
models were constructed using Cloudy ( Fcrl and et al.|[l998l ) , 
version 07.02, and assuming a plane-parallel open geometry 
with density equal to the value adopted fr om the observa- 
tions (see Table [2}. We used a WM-basic l|Pauldrach et al.l 
1200 j ) model stellar atmosphere, with T eff = 39000 K 
and log q = 4.0, value s very similar to those derived by 
ISimon-Dfaz et alj (|2006t l for 6 1 Ori C. The number of ioniz- 
ing photons entering the ionized slab was specified using the 
ionization parameter, U, the ratio of hydrogen-ionizing pho- 
tons to the hydrogen density. We changed the value of this 
parameter till the degrees of ionization (given by + /0 2+ ) 
derived for the models were close to the observed ones. The 
final values used for U imply similar numbers of hydrogen 
ionizing photons in the nebular and shock models, with a 
difference of ~40%. 

We have used the Cloudy "H H region" abundances , 
based on the abund ances derived b y Baldw in et al.1 (|l99ll ). 
iRubin et all (|l99ll ) and lOsterbrock et ah I (|l992h in the 
Orion Nebula, for all elements except iron, which we have 
rescaled in order to reproduce the observed [Fe ill] 4658/H/3 
line ratios. The models also have "Orion" type dust: graphite 
and silicate grains with the Orion size d istribution (deficient 
in small grains, see lBaldwin et al]|l99lf) and an original dust 
to gas mass ratio of 0.0055. For each model, we used the cal- 
culated line intensities to derive the physical conditions and 
the + and 2+ abundances, following a similar procedure 
to the one used to derive the observational results. 
We computed three shock models. Model A has similar char- 
acteristics to the nebular model except for the density, and 
can be used to assess the influence of this parameter on the 
electron temperatures. Model B has the same input param- 
eters as model A, but with the Fe abundance multiplied by 
a factor of 7. Model C has the same input parameters as 
model B but with half the amount of dust. The input model 
parameters are listed in Table [5] Table [§] shows a compari- 
son for observations and models of the physical conditions, 
the degree of ionization given by + /0 2+ , and the [Fe ill] 
4658/H/3 line ratio. We can see that the nebular model re- 
produces well the observational constraints. As for the shock 
models, the increment in the density of model A is enough 
to explain the temperatures found for the shock component, 
but does not reproduce the [Fe III] flux, whereas the higher 
Fe abundance of model B reproduces well this flux. Model 
C illustrates that a reduction in the amount of dust does 
not change significantly the values of the chosen constraints. 
The introduction of grains smaller than the ones considered 



in the Orion size distribution of Cloudy will lead to higher 
temperatures through photoelectric heating, but we do not 
know what grain size distribution would be suitable for the 
shock component. 

The iron abundances in those models that reproduce the 
observed [Fe III] line fluxes are somewhat higher than the 
ones derived from the observations, but they reproduce the 
value of the shock to nebular abundance ratio. This can 
be considered a confirmation of the increment in the iron 
abundance in the shock component, which is most probably 
due to dust destruction. As it has been discussed in H5.6I 
dust destruction could also increase the gaseous abundances 
of other elements like carbon or oxygen and this increment 
will change the amount of cooling and hence the electron 
temperature. We ran a model where the abundances of O 
and C were increased by 14% and found that the derived 
temperatures decrease by about 200 K. 



6 CONCLUSIONS 

We have obtained deep echelle spectrophotometry of 
HH 202-S, the brightest knot of the Herbig-Haro object 
HH 202. Our high-spectral resolution has permitted to sep- 
arate two kinematic components: the nebular component 
-associated with the ambient gas- and the shock one - 
associated with the gas flow. We have detected and mea- 
sured 360 emission lines of which 352 lines were identified. 
We have found a clear disagreement between the individual 
c(H/3) values obtained from different Balmer and Paschen 
lines. We outline a possible solution for this problem based 
on the effects of Ly-continuum pumping and 1-changing col- 
lisions with protons. 

We have analyzed the ionization structure of HH 202- 
S concluding that the dominant excitation mechanism in 
HH 202-S is photoionization. Moreover, the dependence of 
the I(A)sh//(A) no b ratios and the ionization potential, the 
comparison of the spatio-kinematic profiles of the emission 
of different ions, as well as the physical separation estimated 
of 0.14±0.05 pc between HH 202-S and 6 1 On C indicate that 
an ionization front is trapped in HH 202-S due to compres- 
sion of the ambient gas by the shock. 

We have derived a high n e , about 17000 cm" 3 , and similar 
T e for the low and high ionization zones for the shock com- 
ponent, while for the ambient gas we obtain an n e of about 
3000 cm -3 , and a higher T e in the low ionization zone than 
in the high ionization one. We have estimated that the pre- 
shock gas in the inmediate vicinity of HH 202 has a density 
of about 700 cm -3 , indicating that the bulk of the emis- 
sion of the ambient gas comes from the background behind 
HH 202. 

We have derived chemical abundances for several ions and 
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elements from the flux of collisionally excited lines (CELs). 
In particular, we have determined the Ca + and Cr + abun- 
dances for the first time in the Orion Nebula but only for the 
shock component. The abundance of C 2+ , + and 2+ have 
been determined using recombination lines (RLs) for both 
components. The abundance discrepancy factor for 2+ , 
ADF(0 2+ ), is 0.35 dex in the shock component and much 
lower in the ambient gas component. 

Assuming that the ADF and temperature fluctuations are 
related phenomena, we have found a t 2 (0 2+ ) of 0.050 for the 
shock component and 0.016 for the nebular one. The high 
t 2 value of the shock component produces some apparent 
inconsistencies between the total abundances in both com- 
ponents that cast some doubts on the suitability of the t 
paradigm, at least for the shock component. However, the 
fact that the values of the t 2 parameter determined from 
the analysis of the He I line ratios are in complete agree- 
ment with those obtained from the ADF(0 2+ ) supports that 
paradigm. 

Finally, the comparison of the abundance patterns of Fe and 
Ni in the nebular and shock components and the results 
of photoionzation models of both components indicate that 
a partial destruction of dust grains has been produced in 
HH 202-S after the passage of the shock wave. We estimate 
that the percentage of destruction of iron dust particles is 
of the order of 30%-50%. 
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